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ABSTRACT 

We report the discovery of 19 hot, evolved, massive stars near the Galac- 
tic center region (GCR). These objects were selected for spectroscopy owing to 
their detection as strong sources of Paschen-a emission-line excess, following a 
narrowband imaging survey of the central 0?65 x 0?25 (I, b) around Sgr A* with 
the Hubble Space Telescope. Discoveries include six carbon-type (WC) and five 
nitrogen-type (WN) Wolf-Rayet stars, six supergiants, and two B supergiants. 
Two of the supergiants have X-ray counterparts having properties consistent 
with solitary stars and colliding-wind binaries. The infrared photometry of 17 
stars is consistent with the Galactic center distance, but two of them are located 
in the foreground. Several WC stars exhibit a relatively large infrared excess, 
which is possibly thermal emission from hot dust. Most of the stars appear 
scattered throughout the GCR, with no relation to the three known massive 
young clusters; several others lie near the Arches and Quintuplet clusters and 
may have originated within one of these systems. The results of this work bring 
the total sample of Wolf-Rayet stars in the GCR to 92. All sources of strong Pa 
excess have been identified in the area surveyed with HST, which implies that 
the sample of WN stars in this region is near completion, and is dominated by 
late (WNL) subtypes. The current WC sample, although probably not complete, 
is almost exclusively dominated by late (WCL) subtypes. The observed Wolf- 
Rayet subtype distribution in the GCR is a reflection of the intrinsic rarity of 
early subtypes (WNE and WCE) in the inner Galaxy, an effect that is driven by 
metallicity. 
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Introduction 



The Galactic center region (GCR) has the highest star formation rate density in the 
Milky Way, more than orders of magnitude above the average value of the Galactic disk 
( IFiger et al.ll2004j ). This is considered to be a consequence of the large reservoir of molecular 
gas that fills the central half-kiloparsec of the Galaxy, where ~10% of the Galaxy's molecular 
gas is concentrated. With respect to the physical properties that are typical of molecular gas 
further out in the Galactic disk, the GCR medium is significantly hotter, more turbulent, 
and more highly magnetized. It is also subject to mechanical shearing by the strong tidal 
field, and re peated compression from clou d-cloud collisions, stellar winds, and supernova 
shock waves ( iMorris fc Serabyn et al.lll996l ). It has been suggested that the environmental 
conditions in the GCR result in an elevated Jeans mass and an enhanced rate of massive-star 
production. 

At least three extraordinary clusters of massive stars have forme d within the unique 
GCR medium during the pa st several Myr, including the Arches cluster (INagata et al.lll995[) 
and the Quintuplet cluster ( Okuda et al. 1990 ; Nagata et al. 1990 ; Glass et al. 199ol ). both 
of which lie ~30 pc, in projectio n, from the c e nter. Within the inner few parsecs of the 
Galaxy lies the Central cluster (IKrabbe et al.l Il99ll ) , which is gravit ationally bound to 
the Galactic black hole, Sgr A* JSchodel et aD I2OO2I : Icihez et alJ l20oi ). All three of the 
known clusters are massive (M ~ 1O 4 M ), dense, containing h undreds of O star s , and 



dozens of evolved W olf-Rayet stars (WRs) and OB supergiants (IFiger et al.l Il999al 12002 



Paumard et al.l 120061 ) . The collection of three such extraordinary clusters in the GCR is 
evidence favoring the hypothesis that the predominant mode of star formation in the region 
gives rise to dense, massive clusters. However, surveys aimed at identifying massive stars 
throughout the greater GCR have revealed a population of several doze n massive stars that 



have no apparent spatial association with the known s t ellar clusters (|Cotera et al. 



Homeier et al. 



2003 



Mauerhan et al. 



Muno et al.l200a Mikles et al.ll2006l : Mauerhan. Muno. fc Morris 



1999 



2007 



2010all2010bl ). The discovery of these "isolated" massive stars has motivated 



the continued search for additional such cases, with the goal of determining the total extent 
of the inter-cluster population of massive stars and their physical relation to, or independence 
from, starburst clusters like the Arches and Quintuplet. The demography and kinematics of 
these isolated stars will yield valuable insight into the mode of massive star formation in the 
GCR and the dynamical evolution of massive stellar clusters within galactic nuclei. 

The identification of isolated massive stars in the GCR is a difficult task. Heavy extinc- 
tion limits observations of the Galactic center starlight to the infrared part of the electro- 
magnetic spectrum. Infrared measurements sample the Ray leigh- Jeans tail of stellar spectral 
energy distributions where there is a degeneracy of stellar color with effective temperature. 
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This renders early-type and late-type stars photometrically indistinguishable. Narrowband 
photometry aimed at detecting emission-line excesses from hot supergiants can be very effec- 
tive but requires a highly stable point-spread function (PSF) in order to accurately subtract 
the stellar continuum emission from the spectral line emission. This is difficult to achieve in 
seeing limited conditions, especially when extreme stellar confusion is an issue, as it typically 
is when observing the GCR. Therefore, infrared narrowband surveys are best conducted from 
the stable environment of space, where low-background, diffraction-limited observations can 
be executed without the nuisance of PSF variability. Furthermore, space-based narrowband 
observations enable measurement of the Paschen-a (Pa) transition (n=4-3) of H I, which is 
the strongest emission line of this species in the infrared. The A1.87 /im wavelength of this 
transition is also shared with several transitions of He I and He II, which is important in the 
interest of detecting hydrogen-poor WRs that exhibit strong helium emission. 

In this work, we report discoveries of "isolated" emission-line stars in the GCR via 
Pa narrowband imaging from space and subsequent infrared spectroscopy from the ground. 
The new additions add to our recent identification of a new luminous blue variable in the 



GCR ( jMauerhan et al.ll2010bl ). which was the first reported discovery of the spectroscopic 
survey. The successful detection of these stars highlights the significant advantage that space- 
based narrowband searches us ing the Pa line have over groun d-based surveys centered on 



other atomic transitions (e. g., lFigerlll995l ; iHomeier et al.l 120031 ) which have identified only a 



fraction of the isolated emission-line stars. We present here the spectra of new identifications, 
their photometric properties, and discuss the implications that these objects have for massive 
star formation in the GCR. 



Observations 



2.1. Selection of Pa Emission-line Stars 



Our target sample was assembled following the completion of a Pa narrowband imag- 
ing survey of the GCR, performed with the Hubble Space Telescope a nd the Near-Infrared 
Camera and Multi-Object Spectrometer (NICMOS; IWang et al. boioh . The survey, requir- 
ing 144 orbits of observations, covered the central (l,b) ~ 39 x 15 arcmin 2 around Sgr A*, 
which corresponds t o a total p rojected physical area of ~ 90 x 35 pc 2 , assuming a distance 
of 8 kpc to Sgr A* ( )Reidlll993l and references therein). The NIC3 camera was utilized and 
narrowband images were obtained through the F187N (line; A = 1.874 /im; FWHM=0.019 
/im) and F190N (continuum; A = 1.901 /im; FWHM=0.017 /im) filters. The survey resulted 
in the detection of ^150 point-like sources of Pa line excess (i.e., sources with F187N-F190N 
values >5cx above the mean of all point sources in the survey area; Dong et al., in prep). 
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Approximately half of the sources are located outside of the three known stellar clusters, 
and nearly half of those were unidentified prior to this work. 

We have undertaken a near-infrared spectroscopic campaign to characterize the uniden- 
tified sources of Pa line excess. Focusing on the sources having the strongest line excess and 
the brightest near-infrared counterparts (K s < 12.7 mag), we have spectroscopically identi- 
fied 19 new, hot, massive stars in the GCR. So far, our spectroscopic search has had a 100% 
success rate for detecting emission-line stars. Table 1 lists basic data for these sources, includ- 
ing their positions, corresponding JHK S photometry, and narrowband F187N and F190N 
measurements. The sources are named based on their Galactic coordinates, but will hereafter 
be referred to by their numerical entry in Table 1. The positions a nd JHK S photometry 



were taken from the catalog of the SIRIUS survey of the GCR (see iNishiyama et al.l 12006 



and references therein), which generated a catalog of point sources having a 10<r limiting 
K s magnitude of 15.6. The F187N and F190N measurements are from Dong et al. (2011, 
in preparation). We note that the method of extraction f or the narrowband pho tometry 
has been refined since our report on the LBV GO. 120—0.048 (IMauerhan et al.ll2010bl ). where 
we presented o ur narrowband mea surements for this so urce and the kno wn LBVs qF362 
(aka FMM 362; iGeballe et aDbood ) and the Pistol Star jFiger et al.lll998h : so, we take this 
opportunity to include the refined narrowband measurements for all three of these LBVs in 
Table 1. 



2.2. Spectroscopic Observations 



Spectra of four stars were obtained using the SpeX medium- resolution spectrograph on 
the 3 m Infrared Telescope Facility (IRTF) telescope (IRayner et al.l 120031 ) . located on the 
summit of Mauna Kea in Hawaii. SpeX was used in short cross-dispersed mode (SXD), 
utilizing a slit width of ft' 5, providing a spectral resolution of R = X/5X ~ 1200. Flat-field 
images and a wavelength calibration source were provided by obtaining spectra of continuum 
and argon lamp sources. The data were reduced and extracted using the IDL-based software 
packag e Spextool, speciall y designed for the reduction of data obtained with SpeX on the 
IRTF (jCushing et al-lbooJ ). 



The 3 m United Kingdom Infrared Telescope (UKIRT) on Mauna K ea was used to obtain 
spectra of five stars. The UKIRT 1-5 fim Imager Spectrometer (UIST; iRamsay et al.ll2004l ) 
was used in service mode as part of the UKIRT Service Programme. The short-i^ grism and 
4-pixel slit (ft'48) were used, providing a spectral resolution of R ~ 2133 and a wavelength 
range of A2. 01-2. 26 /mi. Flat-fields and a wavelength calibration source were provided by 
obtaining spectra of continuum and argon lamps. The data images were reduced using the 
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Starlink ORACDR pipeline. The spectra were extracted using the IRAF routine APALL. 



Spectra of four stars were obtained using the 4.1 m Anglo- Australian Telescope (AAT) 
on Siding S pring Mountain (M ount Woorat) in New South Wales, Australia. The IRIS2 
instrument (ITinney et al.l 120041 ) provided a spectral resolution of R ~ 2400 in the K band, 
using the 1" slit. Flat-field images and a wavelength calibration source were provided by 
obtaining spectra of continuum and neon lamp sources. The data images were reduced and 
spectra extracted with the same programs used for the UKIRT data. 

Spectra of six stars were obtained at the Southern Observatory for Astrophysical Re- 
search (SOAR), loc ated on Cerro Pac tion in Chile. The Ohio State Infrared Imager Spec- 
trometer (OSIRIS; iDepoy et al.l Il993l ) was used in cross-dispersed mode for five of these 
stars. This mode provides a spectral resolution of R ~ 1200. Wavelength calibration was 
achieved using the background OH emission lines from the sky. For one star, OSIRIS was 
used in the higher resolution longslit mode to obtain a .fT-band spectrum with a resolution 
R ~ 3000. For this observation, wavelength calibration was performed using spectra of the 
internal neon and argon lamps. All of these data were reduced and extracted using the 
standard IRAF routines and APALL. 



For all observations the spectra were acquired in an "ABBA" nodding sequence in order 
to subtract the sky background, and to suppress the contribution of bad pixels. Spectra 
of A0V standard stars were obtained at airmasses similar to those of the science targets in 
order to derive a telluric ab sorption spectrum . The telluric corrections were applied using 
the IDL package xtellcor (IVacca et al.l 120031 ) . which removes model H I absorption lines 
from the A0V standard star before application to the science data. Table 2 summarizes the 
facilities used for the observations of each star. 



3. Spectroscopic Classification 



In the following sections, the spectroscopic criteria we use to classify these stars were 
ad opted from pr e vious s pectroscopic studies of massive sta rs in t he near-infrared, con ducted 
bv iMorris etaD Jl996h . iFiger. McClean k Naiarro Ul997h . and iMartins et al.l (bood ). The 
central wavelengths of all spectral lines referred to in the text are adopted from these refer- 



ences. 
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3.1. Early O Supergiants and Weak-lined WNh Stars 

The .fT-band spectra of stars 6 and 7 are presented in Figure 1. The spectra are domi- 
nated by a complex of blended emission lines near A2. 112-2. 115 /im, which includes contri- 
butions from He I, N ill, C ill and O ill. Blueward of this line complex are emission lines 
of C IV at A2.069 and A2.078 jum. Br7 emission is marginally detected in the spectrum of 
star 6, although it may be a background nebular feature. The spectra of both stars exhibit 
He II absorption at A2.189 fim, although it is very near the noise level. These features are 
consistent with the if-band sp ectra of early O supergiants, specifically in the range of 04-61 



e. g., see iMartins et all 120081 ) . 



The spectra of stars 16, 15, 12, 9, and 10 are presented in Figure 2. The spectra of these 
stars are dominated by emission lines of Br7 and the A2. 112-2. 115 /mi complex of He I, N 
III, C ill and O ill. Stars 16, 15, and 12 exhibit the A2. 112-2. 115 /mi complex in fairly strong 
emission, while 9 and 10, instead, exhibit a He I absorption line accompanied by a relatively 
weak emission component. Weak C IV emission is also apparent in the spectra of stars 16, 
15, and 12. The Br7 lines of all of these stars exhibit an asymmetry on their blue side owing 
to a contribution from He I at A2.1647 /im. An absorption line of He I at A2.058 /im is 
present for four of the stars, while this feature appears in emission for two of them. Weak N 
ill emission appears at A2.247 /im in the spectrum of star 16. Three of the stars also exhibit 
weak He II absorption at A2.189 /im, while this transition is not detected in the spectra stars 
9 and 10. These features are characteristics of OIf+0 supergiants and late-type WN stars 
(WNL), which are commonly referred to as hydrogen-rich WRs (WNh stars). Distinguishin g 



between these two types can be very difficult indeed. According to IMartins et al.l ( 120081 ). 
OIf + stars exhibit Br7 emission that has approximately equivalent strength to the A2.112- 
2.115 /im complex, while WN8-9h stars exhibit a relatively dominant Br7 emission line. 
Based on the dominance of Br7 in the spectra of all stars in Figure 2, we may be inclined to 
classify them all as WN8-9h stars. However, the weakness of the A2. 1 12-2 .115 /xm complex 



in the spectra of stars 9 and 10 is more typical of the 04-6If + stars in IMartins et al.l (120081 ) 
than the WN8-9h stars, which always have a strong A2. 112-2. 115 /mi emission complex. 
Furthermore, WN8-9h always exhibit He II emission or absorption at A2.189 /im, but this 
transition is not detected in the spectra of stars 9 and 10. Thus, we classify stars 16, 15, 
and 12 as WN8-9h stars, and assign the 04-6 If + designation to stars 9 and 10. 



1 Historically, the "+" in the OIf + designation is given to O stars exhibiti ng S IV emissio n (AA4089 and 
4116 A), as well as N in and He II A4686 A emission in their optical spectra ( Walbornlll971 ). 



-7- 



3.2. Other OB Supergiants 

The spectra of stars 13 , 18, 5, and 1 are presented in Figure 3. All of these stars exhibit 
prominent Br7 emission. Stars 13 and 18 also exhibit strong He I emission at A2.058 /im 
with a P Cygni profile, and weak He I absorption at A2.112 /im that may also include a 
weak N ill emission component on the red side of the He I absorption line. The Br7 line of 
star 13 exhibits an asymmetry on its blue side owing to a contribution from He I. Star 1 
exhibits He I at A2.058 /im in emission, although it is relatively weak. Star 5 exhibits the 
He I A2.058 /im feature as well, but it is in weak absorption. The spectra of stars 5 and 
1 exhibit weak emission from the Na I doublet near A2.209 /im, and in the case of star 5, 
weak Fe II emission near A2.09 /im. The presence of these low-ionization emission features 
indicates that stars 1 and 5 are probably of the B spectral class, with subtypes in the range 
of B0-B2 (e.g., see Hanson et al. 1996). Their classification as supergiants is justified in 
Section 4.2. By comparison, stars 13 and 18 must be hotter, as indicated by the presence 
of He I absorption, while the strong P Cygni emission lines of He i at A2.058 /im indicate 
relatively strong winds. Based on these features, we conclude that stars 13 and 18 can be 
classified as P Cygni- type O supergiants." 



3.3. Strong-lined WN Stars 



Figure 4 presents spectra of stars 11 and 17. These stars exhibit broad emission lines, 
indicative of fast, extended stellar winds. For star 11 spectral emission features include Br7, 
He I at A2.058 /im and A2.112 /im, and He II at A2.189 /im. The latter of these features 
is probably the superposition of a He II A2.189 /im P Cygni profile and a blend of He I 
emission at A2.185 /im. Weak emission from N ill is also present at A2.247 /im, as well 
as a weak C Ill/N ill blend at A2.104 /im. The weakness of He II relative to Br7 and He 
I at A2.112 /xm, and the strong P Cygni emission component of He I at A2.058 /im, are 
characteristics of late- type, hydrogen-rich, nitrogen- t ype stars, specifically WN8-9h stars 
( TFiger. McClean fe Naiarro Il997t iMartins et al"1l2008f ). 



The dominant emission feature in the spectrum of star 17 is He II at A2.189 /im, which 
indicates that it is of an early WN (WNE) subtype. Weak but broad He I emission at A2.058 
/im is present with a deep and highly blueshifted P Cygni absorption component. The 
large blueshift of this absorption component relative to the peak of the emission component 
highlights the high wind velocity of star 17 (5A=0.015 /^m =>• t> wind rs 2200 km s _1 ) relative to 
star 11 (<5A=0.006 /im =>■ t> wind rs 900 km s -1 ). Furthermor e, the emission l i nes of star 17 are 
exceptionally broad, having FWHM=155 A. According to lCrowther et al.l ( 20061 ). WN stars 
with FWHM(He II A2.189 /im) > 130 A are classified as broad-lined WN stars (WNb). 
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Star 17 satisfies the broad criterion, while the relative strength of He II relative to Br7 and 
He I is consistent with a WN5 spectral type. Thus, we classify star 17 as a WN5b star. 



3.4. WC Stars 

Figure 5 presents the K-band spectra of star 8, 19, 14, 3, 2, and 4. Each of these 
stars exhibit emission from C ill at A2.115 fim, He I, and He II near A2.058 /xm, C ill, 
and C IV near A2. 07-2. 08 /zm, and weak He II emission near A2.189 /im. These spectral 
characteristics are consistent with those of late, carbon-type (WC) WRs. A subtype di- 
agn ostic is provided by the strength ratio of C IV at A2.08 //m relative to C ill at A2.115 



(IFiger. McClean fc Najarro 1119971 ). which is less than unity for all stars, consistent with 



WC9 subtype. 

4. Photometric Properties 
4.1. Infrared Excess 

Luminous WR and O stars generally exhibit infrared excess fro m the free-free emission 



generated in their ioniz ed winds, and also from strong line emission (IHadfield et al. 112007a 



Mauerhan et al.ll2009bl ). Thermal continuum excess from hot dust emission is also not un- 
common to late-type WC (WCL) stars. Color-color diagrams are useful in assessing the 
infrared excess of WRs and the likely causes for it. Figure 6 shows a color-color diagram of 
H — K s versus J — H for the stars in our sample that were detected in all three near-infrared 
bands. The interstellar reddening vector is illustrated by the linear trend traced by random 
field stars that lie toward the GCR, most of which are presumably late-type giants. Stars 1 
and 7 lie near the lower left of Figure 6, well separated from most of the other massive stars. 
The fact that they suffer relatively low extinction implies that they are in the foreground 
toward the GCR. They also do not appear to exhibit evidence for significant, intrinsic excess. 

While the OB and WN stars in the GCR have color excesses typical for their class, 
several WC stars in Figure 6 appear to have color excesses that are significantly higher than 
those of the O and WN stars and 1 WC star (note that only four of the WC stars from our 
sample are included in Figure 6, since two others have not been detected in the J band). 
The colors of star 19 are consistent with the excess being from free-free emission alone, and 
this might also be true for star 14. However, the relatively strong excess emission from 
stars 2 and 4 (the two WC stars that are furthest from the reddening trend in Figure 6) 
indicates that there is an additional excess component, which is likely the result of thermal 
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emission from hot dust. Late-type WC stars (WC8-9) commonly exhibit thermal emission 
from hot dust, and it is generally accepted that the production of dust from these stars is a 
consequence of binary evolution, whereby the dust is created via the collision of hydrogen - 
rich and carbon-rich winds from the respective OB and WC stars (jWilliams et al.l 120051 ). 
Several likely examples of this phenomenon are alre ady known in the GCR , and they are 
typically detectable as hard X-ray sources (e.g., see iMauerhan et al.l l2010al and references 
therein). Since sta rs 2 and 4 were not detected in X-ray surveys of their respective fields 
( jMuno et al.l 120091 ). there is currently no evidence that they are members of colliding-wind 
binary systems. Nonetheless, the relatively strong excess emission from stars 2 and 4 is 
consistent with thermal dust emission. Thus, we tentatively classify stars 2 and 4 and 
WC9?d stars. 



4.2. Extinction and Luminosity 



We derived the extinction and absolute photometry for each source by comparing the 
observed J — K s and H — K s colors with the intrinsic colors of stars of the same spectra l 
types, which were adopted from t he literature (iMartins &z Plezll2006l ; ICrowther et al.l 120061 ) . 
We used the extinction relation of iNishiyama et al. for stars near the Galactic center, 

which is given by 

A k =1AA±0ME h _ Kb 



and 



A Ks = 0.494 ± 0.006£j_^ s 



where E(H — K s ) = (H — K s ) ^, s — {H — K s ) , and similarly for the J — K s colors. The 
two values were then averaged to obtain a final extinction estimate, Ak s - The results are 
listed in Table 3. The values are similar to other known massive stars in the GCR (e.g., see 



Mauerhan et al.ll2010al ). As expected from their positions in Figure 6, stars 7 (04-61) and 1 
(B0I-B2I) have relatively low Ak s values of 1.16 and 0.95 mag, respectively, which indicates 
that they lie in the foreground. 



Assuming a distance of 8 kpc for the Galactic center (lReidlll993l ). we derived M Ks from 
the extinction-corrected values. To calculate the bolometric luminosities of the stars, bolo- 
metric corrections for the appropriate spectral types were adopted from the same literature 
sources used to obtain the intrinsic colors, except for the WN8-9h stars, whose bolomet- 
ric corrections were derived from the Mr and luminosity values of the WN8-9h stars in 
Martins et al.l (120081 ) . These results are also listed in Table 3, and they are more or less 
consistent with luminosities that are typical for these spectral types. The luminosity values 
derived for the other O and B stars are also consistent with supergiant luminosity class. 
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We note that the luminosity values obtained for WC9?d stars should be met with caution, 
since thermal dust emission can have a significant effect on the observed infrared colors, 
leading to an overestimate of the extinction toward the source. Nonetheless, the luminosity 
values that we obtain are consistent with those of other WC9 stars. For the foreground 0-41 
star 7, we assumed that it shares the same value of Mk as star 6, since they appear to be 
identical in spectral type. Thus, we derive a distance m odulus of 12.79 mag (3 .6 kpc) to 
star 7, which places it in the Norma arm of the Galaxy (IChurchwell et al.l 120091 ). For the 
B0I-B2I star 1 we have no information on its luminosity class, so we are unable to derive 
its foreground distance. However, its extinction value is very comparable to star 7; so, we 
could infer that it also lies at roughly the same distance, also in the Norma arm. If so, then 
the associated distance modulus of 12.79 mag indicates that star 1 has Mx a = —6.71, which 
implies supergiant luminosity class. 



5. X-ray Detections 



The 04-61 stars 6 and 7 have X-ray counterparts in the master catalog of lMuno et al.l ( 120091 ). 
which contains 9017 X-ray point sources detected in all combined observations of the GCR 
with the Chandra X-ray Observatory. The X-ray photometric data for t hese two sources i s 
presented in Table 4. The details of the extraction process are described in lMuno et al.l ( 120091 ). 
C andidate near-inf r ared co unterparts to the population of GCR X-ray sources were identified 
in iMauerhan et al.l ( 12009ai ); the 04-61 star 6 is source 432 from Table 3 of those authors, and 
is associated with the Chandra source CXOGC J174531. 4— 285716. On the other hand, the 
foreground 04- 61 star 7 that is a s sociate d with Chandra source CXOGC J174537.9— 290134 
is not present in IMauerhan et al.l (j2009al ) , since it was flagged as having erroneous photome- 
try in the SIRIUS near-infrared catalog that was used for that work; so, for this star we relied 
on photometry from the Two Micron All Sky Survey (2MASS; Cutri et al. 2003). Table 4 
includes X-ray source astrometry, total on-source integration time, photon number counts 
in the hard (2.0-8.0 keV) and soft bands (0.5-2.0 keV), broadband photon flux, hardness 
ratios for soft and hard energy bands, and the average energy of all detected photons. 

In order to constrain the distance and nature of the X-ray sources, it is useful to examine 
the hardness ratio, defined as HR=(/i — s)/(h + s), where h and s are the fluxes in the hard 
and soft energy bands, respectively. The soft color, HRO, is defined by letting h and s be 
the fluxes in the respective 0.5-2.0 keV and 2.0-3.3 keV energy bands. The hard color, 
HR2, is defined by letting h and s be the fluxes in the 4.7-8.0 keV and 3.3-4.7 keV energy 
bands, respectively. Galactic center sources, which typically lie behind a hydrogen absorption 
column of iV H > 4 x 10 22 cm -2 , will have most of their soft X-ray photons absorbed by the 
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interstellar medium (ISM), and will, thus, exhibit HRO > —0.175 ( IMauerhan et al.ll2009al ). 
Star 6, with HR0=0.56, definitely satisfies the criterion for Galactic center location, while 
star 7, with HR0=— 0.17 barely makes the cut. Compared with HRO, the more energetic 
photons used to calculate HR2 suffer significantly less absorption from intervening gas and 
dust, so they are more useful for constraining the intrinsic X-ray hardness. We compared 
the HR2 colors of these sources with th ose of the several dozen X-ray emitting massive 
stars in the GCR that were presented in IMauerhan et al.l (l2010al ; their Figure 11). Star 6 
is among the faintest X-ray sources that has a confirmed infrared counterpart. It was not 
significantly detected at energies above 4.7 keV, so its HR2 value is —1.0. This contrasts 
with most of several dozen known, massive stellar X-ray sources in the GCR, many of 
which have significant detections at energies up to 8.0 keV. This feature of the general 
population is presumed to be the result of emission from a relatively hard thermal plasma 
with kT > 2 keV. By contrast, the practically null detection of X-rays above 4.7 keV energy 
from star 6 implies that it is a relatively faint and soft X-ray source (kT < 1 keV, if it 
is thermal). Unfortunately, there are not enough counts to effectively model the X-ray 
spectrum. But if we assume a thermal plasma and guess a reasonable temperature, we can 
at least get an order-of-magnitude estimate of the X-ray luminosity. To do this, we used the 
derived extinction values in Table 3 to calcu late the corresponding h ydrogen column density 
(iV H ), using the Ay/Ax ~ 16.1 ratio from iNishiyama et al.l (120081) . assuming a ne gligible 
difference between Ak and Ak s , and the Av/N H ratio of iPredehl &: Schmittl (119951 ). which 
gives N-n/Ay = 1.8 x 10 21 cm -2 mag -1 . We obtained iV H = 5.9 x 10 22 cm -2 for star 6, which is 
consistent with the canonical average value of iV H = 6 x 10 22 cm -2 that is typical for Galactic 
center X-ray sources. Using the observed 0.5-8.0 keV energy flux of Fx = 7.3 x 10~ 16 erg s _1 
cm -2 for this source, we employed the PIMMS0 simulator to calculate the absorption-corrected 
X-ray flux. We assumed an intrinsic thermal plasma energy value range of kT = 0.6-1 keV, 
which is appropriate; after all, the faintness of this source and its la ck of high energy emissio n 
(low HR2 value), compared with the other GCR X-ray sources of IMauerhan et al.l (j2010al ). 
imply that the plasma temperature is unlikely to be any hotter than this. PIMMS predicts an 
unabsorbed X-ray flux in the range of (1.7-9.2) x 10 -14 erg s _1 cm -2 for star 6. Accounting 
for the distance of 8 kpc, this implies an X-ray luminosity in the range of (1.3-7.0) x 10 32 erg 
s _1 in the 0.5-8.0 keV energy band. Thus, assuming kT is between 0.6 and 1 keV, the X-ray 
to bolometric luminosity ratio of star 6 is in the range of log (L x /L hol ) = —7.2 to —6.4. 

The other 04-61 star, 7, has a relatively large X-ray flux in the soft energy range (0.5- 
2.0 keV) compared with star 6 and the other massive stellar X-ray sources in the GCR. It 
is also relatively bright overall. This is not surprising, since both the photometry and the 



2 http://heasarc. gsfc.nasa.gov/Tools/w3pimms. html 
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extinction value that we calculated for the stellar counterpart in Table 3 imply a foreground 
distance of 3.6 kpc. So, if the X-ray-detected 04-61 stars 6 and 7 have the same intrinsic 
nature, then we should expect the foreground star 7 to appear brighter. Fortunately, star 
7 has enough counts to warrant a fitting of its X-ray spectrum, which can better constrain 
the physical p arameters of the X-ray emission. To do this, we used the program XSPEC 



( lArnaudlll996[ ). First, the X-ray spectral data were binned to have at least 30 counts in each 



energy bin. Then, we performed least squares fits to the data by trying out two different 
emission models, including a power-law model and a thermal plasma model. The power-law 
fit required an unrealistically steep photon index of T « 5, so we do not regard this as a 
viable model for star 7. However, the thermal plasma model fit the data reasonably well. The 
parameters fit by the model were iV H , kT, and a normalization factor that is proportional 
to the emission measure. For the iterative fitting procedure, we used initial guess values of 
N H = 3.4 x 10 22 cm -2 , derived from the infrared photometry, and kT = 1 keV. After 30 
fitting iterations, XSPEC returned values of N n = (3.0 ±0.4) x 10 22 cm" 2 , kT = 0.9 ±0.2 keV, 
and an energy flux of Fx = 4.2 x 10~ 15 erg s _1 cm~ 2 . The reduced \ 2 °f the fit was 0.7 for 
13 degrees of freedom. Using these parameters we calculated an absorbtion-corrected X-ray 
flux of F x ™ bs = (5.8 ± 1.3) x 10~ 14 erg s" 1 cm -2 . Using the derived distance of 3.6 kpc, we 
calculated an X-ray luminosity of Lx = 1.3 x 10 32 erg s _1 for star 7, which implies an X-ray 
to bolometric luminosity ratio of log (Lx/L bol ) = —7.2, 

The estimated plasma temperatures and X-ray luminosities of stars 6 and 7 are con- 
sistent with the values that have been measured for the GCR 04-61 sta rs in the sample of 



Mauerhan et al.l feOlOaf ). and the 04-61 stars in the Cygnus OB2 sample of lAlbacete Colombo et al.l ( 120071 ) 



O stars typically ex hibit log (Lx/L hn] ) —7, and this is true for both single stars and those 



in massive binaries ( lOskinoval 120051 ) . This relation is a result of the connection between the 
strength of radiatively driven winds and the temperature of the resulting shocks, which may 
occur inside the stellar wind via turbulence-induced micro shocks, or from wind collision 
if the O star is in a massive binary. However, the plasma temperatures reached via the 
inter-wind micro-shock process typically have kT < 1 keV, while shocks in colliding-wind 
binaries generate higher temperatures (kT > 1-2 keV) and produce harder X-ray spectra. 
Our derived value of kT = 0.9 ± 0.2 keV is insufficient to effectively discriminate between 
these two mechanisms. Although the average photon energy of both sources is (E) = 2.5-3.0 
keV, these photons could simply be the tail of a much softer intrinsic spectrum whose lower 
energy photons were absorbed by the ISM. Alternatively, our assumption of a single tem- 
perature plasma for star 7 might be too simplistic, and the highest energy photons detected 
from this source could actually be from an intrinsic hard component with kT > 1-2 keV. But 
without more counts in the softer energy bands, we cannot reliably fit a multi-temperature 
plasma model to the data, which is generally the case for high X-ray flux O stars in Cygnus. 
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In conclusion, the X-ray properties of stars 6 and 7 are consistent with thermal emission 
from either single O stars or those in colliding-wind binaries, although we currently have no 
good reason to invoke binarity for either source. 



6. Discussion 

6.1. How Complete is the WR Sample in the Central ^90 Parsecs of the 

Galaxy? 

There are currently 92 identified WRs in the GCR, including all known WNE (WN4-6), 
WNL (WN7-9), WCE (WC4-7) and WCL (WC8-9) stars identified in the field (including 
those identified outside of the Pa survey area), and those identified in the Arches, Quintuplet, 
and Central clusters. Their subtype statistics are presented in Table 5. The most noteworthy 
characteristic of the WR subtype distribution is the high late-type to early-type ratio for 
both WC and WN stars; out of the 55 WN stars, only 5 are WNE types; and out of 36 WC 
stars, only 1 is a WCE type. WC stars comprise 40% of the total WR population (47% if 
we exclude the Arches cluster, which is too young to have produced WC stars). 

Figure 7 shows a plot of Pa line strength (defined here as the F187N/F190N flux ratio) 
versus iC^-band magnitude for known WRs in the GCR and the remaining unidentified Pa- 
excess candidates in the survey area (jWang et al. boiol : Dong et al. 2011, in preparation). 



The known sources in Figure 7 include emission-line stars from the GCR field, and from 
the Arches, Quintuplet, and Central Parsec clusters. The fC,-band values of unidentified 
sources were obtained by cross-correlating the positions of the Pa-excess candidates with 
sources from the SIRIUS and 2MASS catalogs (we make a note of caution that a few of 
the unidentified Pa sources having the faintest K s matches might be spurious associations 
in this crowded region, and that SIRIUS magnitudes may have significant systematic errors 
because of source confusion, especially near the Central Parsec Cluster). The data in Figure 
7 suggest that all but a few of the remaining unidentified Pa candidates occupy a brightness 
range of K s and F187N/F190N that is systematically offset from the range occupied by 
confirmed WRs. Thus, an important question naturally arises as to the completeness of the 
WR sample at the inner (I, b) ~ 90 x 35 pc 2 around Sgr A* (the survey area). 

WNL stars, which almost completely dominate the current WN sample, have F187N/F190N 
ratios between 1.2 and 2.9, and have bright counterparts with K s pa 7.5-12 mag. WNL types 



are expected to generally have K < 12 mag when observed in the GCR (e.g., see [Figerj [1995 
their Figure 25). Thus, we suspect that only a few WNL stars, at most, could remain 
unidentified in the survey area. This also appears to be the case for WNE types; those 
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confirmed thus far, all of which have subtypes in the WN5-6 range, have several of the high- 
est F187N/F190N ratios in the entire sample of WRs. None of the unidentified Pa-excess 
candidates have F187N/F190N ratios this large. Earlier WNE stars of subtype WN3-4 have 
not yet been found in the GCR. Although they are intrinsically the faintest of the WN stars 
(Mk ~ —3.1 mag; ICrowther et al.ll2006l ). having expected K-baiad magnitudes as faint as 
f=sl4 in the GCR, WN3-4 stars will exhibit Pa excess similar to the WN5-6 stars. Therefore, 
they should still be easily detectable as relatively strong sources of Po excesfl modulo the 
presence of a bright, diffuse background. Thus, although Figure 7 indicates that several WN 
stars could remain unidentified, we suspect the WN sample in the survey area is near comple- 
tion nonetheless. This implies that WNL types heavily dominate the GCR WN population, 
and that WNE subtypes are intrinsically rare in the region. 

The situation is less clear for WC stars. Their F187N/F190N ratios span a range 
of values similar to the distribution of WNL types, but the faint end of their expected 
brightness distributio n in the GCR extends down to K « 15 mag for the earliest WCE 
subtypes (iFigerl 119951 ). This has been c onfirmed by photome try of the WC5-6 star IRS 3E 
of the Central Parsec cluster (E58 from lPaumard et al.ll2006l ). which has K = 15.0 mag and 
is the only WCE star identified in the GCR to date. Unfortunately, the F187N/F190N ratio 
for this source could not be measured because of source confusion and the strong, diffuse 
background emission from Sgr A West. Even with relatively large intrinsic F187N/F190N 
ratios, however, sources as faint as IRS 3E will have large photometric uncertainties that 
could render them weak or undetectable sources of Pa excess in the HST /NICMOS survey 
(Dong et al. 2011, in preparation). Thus, it is difficult to ascertain the detectability of a 
hypothetical population of WCE types in the GCR, and we cannot reliably estimate the 
potential contribution of such stars to the sample of unidentified sources in Figure 7. The 
same is true for WCL types, which can also be very faint (K < 14 mag) and whose line 
strengths can be heavily diluted in cases where there is signi ficant flux contributio n from a 
host dust continuum. Extremely dusty WCL stars (DWCL; [Williams et al.lll987l ). like the 
Quintuplet proper members, whose emission lines have been completely diluted by a thermal 
continuum or obscuring dust screen, are not detectable Pa excess sources at all; so, some of 
these stellar types could remain unidentified. However, DWCL stars are very bright infrared 
sources (K a ~ 6-7 mag) and typically have X -ray counterparts . Based on these properties, 
the statistics of X-ray /infrared matches from iMauerhan et all (j2009al ). combined with the 



3 We have confirmed this statement by examining the HK spectra of WN3-6 and WC4-6 stars obtained 
from the ground (not shown here). The spectra of the WNE stars exhibit iarge flux at A1.87 /im that can 
penetrate the severe atmospheric absorption at this wavelength, which implies that such stars will generally 
have high F187N/F190N ratios. 
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spectroscopic survey results of iMauerhan et al.l (l2010al ). imply that most, if not all, DWCL 
stars in the survey area should have been found already. Still, less extreme WC9d stars with 
diluted emission lines that are not detectable X-ray sources might remain elusive. 

The near completeness of the WN stars in the survey area implies that apparent domi- 
nance of WNL stars over WNE stars is an intrinsic property of the GCR. Although the sample 
of WC stars is probably not complete, we also suspect that WCL stars are intrinsically dom- 
inant over WCE types. The reasons for this are as follows. It has been well established 
that WR subtype distributions are correlated with environmental met allicity; this h as been 
established for the Milky Way and other galaxies of the Local Group (iMasseyl 119961 ). WNE 
and WCE stars appear to be rare in relatively high-metallicity environments (Z > Z ), such 
as the inner Galaxy, but are common in low-metallicity environments, such as the outer 
Galaxy [Z > Z & ), and the Large Magellanic Cloud (Z <C Z Q ). This observed trend has been 



interpreted as a result of the metallicity dependence of WR winds (see ICrowtherl 120071 and 
references therein for the following discussion). For WN stars, the higher wind density that 
results from increased metallicity leads to more efficient atomic recombination from high 
ionization stages (e.g., N 5+ ) to lower ions (e.g., N 3+ ), which effectively increases the trend 
toward later spectral types. Alternatively, decreased recombination efficiency and reduced 
nitrogen content at lower metallicity results in a trend toward earlier WN subtypes. For WC 
stars in high-metallicity environments, higher wind densities lead to an increased emission 
strength of C ill with respect to C IV. Since WCL and WCE stars are classified as such 
based upon the relative strengths of C ill and C IV emission lines, WC stars in environments 
with Z > Z & will tend to manifest as WCL subtypes ( C ill stronger than C iv). Although 
this effect was originally suggested by ICrowther et al.l (120021 ) to be an important factor for 
the optical classification of WC stars (specifically for C ill at A5696 A), the effect of elevated 
C ill relative to C IV for increasing wind density is likely to be mirrored in the infrared 
spectrum of WC stars as well. Thus, the dominance of WNL and WCL stars over WNE 
and WCE stars in the GCR is a natural conseque nce of the fact that the metallicity of the 



flCunha et all 120071 : iNaiarro et al.l 120091 ) 



region is at least solar fe.g.. seelDavies et al.ll2009l and references therein), if not supersolar 



6.2. Rogue Massive Stars from the Arches and Quintuplet Clusters 

Figures 8 and 9 show the distribution of "isolated" massive stars very near the Arches 
and Quintuplet clusters. For the Arches, there are three WNh stars that lie within ^2 
pc, in projection, of the cluster. The magnitudes and spectral types of these outer WNh 
stars are very similar to the brightest Arches members, some of which have estimated ini- 
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tial masses above 100 M Q (IFiger et al.ll2002l ). Based on the spectral similarities and the 
close spatial proximity, we suggest that the three massive stars surrounding the Arches 
cluster are coeval with the system. The main uncertainty is with regard to their current 
kinematical state. One might expect that mass segregation via dynamical friction should 
concentrate massive stars like these toward the cluster center. However, higher order pro- 
cesses occurring at the cluster center might result in an alternative effect. Indeed, the Arche s 
cluster is the most dense system known in the Galaxy (3 x 10 5 M & pc -3 ; IFiger et al.lll999bl ). 
As such, the central region of this cluster must be the site of frequent gravitational en- 
counters between massive stars and binaries, including binary captures, disruptions, and 
exchanges. Through such interactions, one might expect a few massive stars to be grav- 
itationally "kicked" fro m time to time beyond the tidal radius, which is ~1 pc for the 

Arches and Quintuplet ( Portegies Zwart et al. 12002 ). T his co uld be the case for stars 11 

and 12 of this work, and GO. 10+0.02 from lCotera et al. I ( 119991 ). all of which lie well beyond 
the Arches cluster's tid al radius. Given t hat the internal velocity dispersion of the Arches 
cluster is ~10 km s _1 ( jFiger et al.l 120021 ). a massive star would not require an unusually 
large gravitational "kick" to escape the cluster. Alternatively, one or more of these stars 
could have formed at the same time and from the same cloud as the main cluster but never 
became bound to the system. If this is the case then the stars could simply be follow- 
ing a path that is in tandem with that of the cluster, without orbiting the cluster center. 
These p ossibilities may also be appli c able t o the massive stars surrounding the Quintuplet 
cluster. iMauerhan. Muno. fe Morris I ( 120071 ) suggested that the 06If + X-ray source CXOGC 
J174617.0— 285131, a suspected colliding-wind binary, may have ori ginated in the Quintuple t, 
while the same might be true for the nearby LBV GO. 120—0.048 ( IMauerhan et al.ll2010bf ). 

Observationally distinguishing between the various possible origins of the stars surround- 
ing the Arches and Quintuplet clusters would be very difficult using radial velocities alone, 
since both ejected stars and comoving stars that were never bound to the clusters would ex- 
hibit similar radial velocities, perhaps differing on the order of only 10 km s _1 . Proper motion 
measurements, however, would allow for the derivation of a three-dimensional space velocity 
vector, which would better elucidate the association between these stars and the clusters, and 
could provide evidence for ejection, specifically. Such measurement s are within the cap abili- 
ties of 10 m class telescopes utilizing adaptive optics. For example, IStolte et al.l ( 120081 ) used 
proper motion measurements over a 4.3 year baseline to derive a three-dimensional space 
velocity of 232 ± 30 km s _1 for the Arches cluster. As coeval siblings of the Arches cluster 
stars, we might expect stars 11, 12, and GO. 10+0.02 to have a similar velocity, whether they 
were ejected from the cluster or are just comoving with it. Multiple epoch measurements 
with longer time baselines, however, could discriminate between ejection and comoving sce- 
narios by achieving the required accuracy to obtain a velocity difference vector between the 
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individual star vectors and the average cluster vector. If the stars were ejected from the 
cluster, we would expect these difference vectors to point away from the cluster. 



6.3. Are there Unidentified Stellar Clusters in the GCR? 



The continuum-sub tracted Pa mosaic of the inner (I, b) ~ 90 x 35 pc 2 of the GCR 
from IWang et al. 1(120 101 ) is shown in Figure 10, where the positions of the known clusters 
and the "isolated" massive stars are illustrated. Aside from the several stars potentially 
associated with the Arches and Quintuplet clusters, some of the other new discoveries appear 
to reside in regions that were heretofore unidentified as potential associations. There is a 
noteworthy association of 3 WC9 stars, a WN7, and an OB supergiant at (I, b) m (359?68 — 
— 359?76, — 0?05). Several promising unidentified candidates having K s rj 12.4 and ~ 12.7 
mag and F187N/F190N ratios consistent with them being WC or OB stars also lie in this 
region. This area is far from any of the known clusters, so, these stars might be the brightest 
representative members of a bona fide association. There is a more highly concentrated group 
containing a WC9d, an Ofpe/WN9, and a B2 la star at (lb) sa (359?86, 0? 00), originally 
identified as a potentially coeval system by lMauerhan. Muno. &: Morris 1 ( 120071 ) . In addition, 
there appears to be a noteworthy collection of 2 WC9 stars and an O supergiant near 
(/, b) « (0?07, -0?06), and a significant concentration of WN stars near (I, b) « (359?98, 0?00) 
and (I, b) ~ (0?05, 0?0 2), including the central stars of the radio H II regions HI, H2, H5, and 
H8 rtzhao et al.lll993h . 



Could these associations of massive stars represent the "tip of the iceberg" of stellar 
clusters that are less compact than the Arches and Quintuplet? There is no measured stellar 



density enhancement in the vicinity of isolated massive stars in the GCR (IWang et al. 112010 
Dong et al. 2011, in preparation). But because of the strong tidal field in the GCR, 
stellar clusters li ke the Arches and Quintuplet will dissolve on a timescale of ~10 Myr 
( IKim et all 119991 ). while the time it takes for less extreme stellar cluster s to fa de into the 
background may be shorter than this. According to iPortegies Zwart et al.l ( 120021 ) the cluster 
dissolution timescale is a sensitive function of galactocentric radius, initial cluster mass, 
and initial density, such that less compact stellar clusters could become invisible as stellar 
density enhancements on a shorter timescale than the Arches or Quintuplet ages, or could 
even persist as undetectable entities for their entire lives. In such cases, an Earth-based 
observer undertaking a narrowband imaging survey of the GCR might see a small group of 
evolved emission-line stars with no obvious surrounding stellar cluster. This could be the 
case for the associations that we have noted above. 



As an alternative explanation for the "isolated" massive stars, there may be two modes of 
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massive star formation operating in the GCR, one that produces dense starburst systems like 
the Arches and Quintuplet, and one that produces massive stars in small, sparse associations 
or in total isolation. Currently, we cannot discriminate between cluster and isolated star 
formation origins; but continued surveys for coeval massive stars in the vicinity of the new 
emission-line stars could provide inroads to solving this problem. 

The project is partly supported by NASA through a grant HST-GO-11120 from the 
Space Telescope Science Institute, which is operated by the Association of Universities for 
Research in Astronomy, Inc., under NASA contract NAS 5-26555. 
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Note. — Positions and JHK S photometry are from IRSF SIRIUS observations of the Galactic center (Nishiya ma et al. 2006), except 
for star 1 and 7, which correspond to 2MASS J17445945-2905258 and 2MASSJ17453798-2901345, respectively iCutri et alj|2003h . The 
HST/NICMOS narrowband F187N and F190N measurements are from Dong et al. (in prep.). 

a This source is the LBV GO. 120-0.048 from iMauerhan et al.l feOlObT) . The JHK a values for this source are from 2MASS. The 
narrowband measurements reported here are refined values that were extracted from the same data, owing to a recent improvement 
in the photometry ext raction method of Dong et al. (2011, in preparation). In addition, the F190N and F190N/F187N measurements 
that were presented by Maucrhan et al. | l|2010bl) for the Pistol Star and qF362 should be updated to 988 ± 18 mjy and 1.21 ± 0.03 for the 
Pistol Star, and 740 ± 13 mjy and 1.02 ± 0.03 for qF362. 
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Table 2. Spectroscopic Observations 
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Telescope /Instrument 
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Table 3. Extinction and Absolute Photometry for Selected Stars 



Star Spectral K s (J — K 8 ) (H — K S ) Q A J K ~ Ks A^,~ Ks A Ks M Ks BC K a L bol 





Type 


(mag) 


(mag) 


(mag) 


(mag) 


(mag) 


(mag) 


(mag) 


(mag) 


(log Lq) 
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B0I-B2I 


7.03 


-0.21 


-0.10 


0.94 


0.96 


0.95 b 


-6.71 b 


-3.7 


6.06 b 


2 


WC9?d 


11.11 


0.23 


0.26 


2.91 


2.86 


2.89 


-6.28 


-3.6 


5.85 b 


3 


WC9 


12.70 


0.23 


0.26 




4.06 


4.06 


-5.86 


-3.6 


5.68 
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WC9?d 


11.14 


0.23 


0.26 


2.21 


2.15 


2.18 


-5.54 


-3.6 


5.56 c 
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B0I-B2I 


11.36 


-0.21 


-0.10 




3.16 


3.16 


-6.30 


-3.7 


5.90 
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04-61 


11.37 


-0.21 


-0.10 


1.96 


2.05 


2.01 


-5.14 


-4.3 


5.68 
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04-61 


8.90 


-0.21 


-0.10 


1.17 


1.14 


1.16 b 


-5.14 b 


-4.3 


5.68 b 
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WC9 


12.02 


0.23 


0.26 




3.12 


3.12 


-5.60 


-3.6 


5.58 
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04-6If+ 
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-0.21 


-0.10 


2.54 


2.67 


2.61 


-7.25 


-4.3 


6.52 


10 


04-6If+ 


10.11 


-0.21 


-0.10 


2.37 


2.40 


2.39 


-6.78 


-4.3 


6.33 
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WN8-9h 


11.14 


0.13 


0.11 




3.41 


3.41 


-6.77 


-4.2 


6.29 


12 


WN8-9h 


11.00 


0.13 


0.11 


2.09 


2.05 


2.07 


-5.57 


-4.2 


5.81 


13 


P Cyg-type OI 


9.12 


-0.21 


-0.10 


2.24 


2.41 


2.33 


-7.71 


-3.7 


6.46 


14 


WC9 


11.61 


0.23 


0.26 


2.41 


2.14 


2.28 


-5.17 


-3.6 


5.41 


15 


WN8-9h 


10.79 


0.13 


0.11 


2.01 


1.98 


2.00 


-5.71 


-4.2 


5.86 


16 


WN8-9h 


10.46 


0.13 


0.11 


2.17 


2.24 


2.15 


-6.19 


-4.2 


6.06 


17 


WN5b 


11.43 


0.36 


0.26 


2.40 


2.09 


2.24 


-5.31 


-4.4 


5.54 


18 


P Cyg-type OI 


9.55 


-0.21 


-0.10 


3.37 


3.35 


3.36 


-8.31 


-3.7 


6.70 


19 


WC9 


11.34 


0.23 


0.26 


2.10 


1.92 


2.01 


-5.17 


-3.6 


5.41 



a Bolometric corrections were derived fromlMartins et al.1 (120081 1 for WN8-9h stars; from lCrowther et alj 1120061 1 
for other WRs; and from lMartins fc Pled (120061') for the OB stars. 

b The relatively low value of extinction for the 04— 61 star 7 and B0I-B2I star 1, compared to the other stars in 
this table, implies that they are located in the foreground. We estimated the distance to 7 (04— 61) by adopting 
Mk s and Lbol from the similar 04— 61 star 6 . The resulting distance modulus is 12.79 mag (3.6 kpc), which implies 
that 7 lies in the Norma arm of the Galaxy l lChurchwell et al]|2009h . Since 1 has an extinction that is comparable 
to 7, one might also assume that these stars have comparable distances. If so, this implies 1 has Mk 3 = —6.71 
mag, which means it is a supergiant. 

c These values may be erroneously high. The color-color diagram in Figure 6 implies that these stars may emit 
thermal excess from hot dust, in which case the extinction cannot be reliably derived using the near-infrared colors. 
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Table 4. Chandra X-ray Data for G359.973-0.008 and G359.925-0.049 



Value G359.973-0.008 (Star 6) G359.925-0.049 (Star 7) 

Spectral type 04-61 04-61 

CXOGC J 174531.4-285716 174537.9-290134 

R. A. (J2000) 266.38108 266.40829 

Decl. (J2000) -28.95466 -29.02625 

a x (arcsec.) 0.5 0.3 

C n s c f (0.5-2.0 keV) 10-8±?;1 140.6jjj;!| 

C n T d (2.0-8.0 keV) 58.81™ 150.8^ 

F tot (s" 1 cm- 2 ) 2.44 x 10~ 7 1.37 x 10~ 6 

HRO 0.57±g;|f -0.17i^o 

HR2 -1.00 -0.07^26 

(£)/photon (keV) 3.0 2.5 

kT (keV) 0.6-1.0 0.9±0.2 

iV H (10 22 cm- 2 ) 5.9 3.0 ± 0.4 

Norm, factor • • ■ 3.2 ± 2.3 x 10~ 5 

F un,bs ( erg s -i cm -2) (1.7-9.2) x 10~ 14 5.8 ± 1.3 x 10~ 14 

log L x (erg s" 1 ) (1.3-7.0) x 10 32 1.3 x 10 32 

log (L x /L hol ) -7.2 to -6.4 -7.2 



Note. The following X-ray data were taken from the catalog of 

Muno et al.l (120091 ): X-ray counts in the soft (C n s c ° ft ) and hard (C n ^ rd ) energy 
bands, hardness ratios HRO and HR2, photon flux F tot , and average photon en- 
ergy (E). All other values were derived here. For G359. 973— 0.008 (Star 6), the 
values of kT, N H , and F x ™ hs were derived from the infrared and X-ray photome- 
try. For G359. 925— 0.049 (Star 7), these sam e values were d erived via a thermal 
plasma model fit, using the program XSPEC ( lArnaudlll996l ). which also yielded 
the model normalization factor. Uncertainties are present when available. 



Table 5. Distribution of WR Subtypes in the GCR 



Group 


WNE 


WNL 


WCE 


WCL 


WR tot 


WC/WN 


WCL/WC 


WC/WR tot 


References 


GCR "Isolated" 


3 


15 





10 


28 


0.56 


1.00 


0.36 


1, 2, 3, 4, 5, 6, 12 


Arches 





13 








13 


n/a 


0.00 


0.00 


7, 11 


Quintuplet 


1 


5 





11 


20 


2.33 


1.00 


0.70 


2, 4, 8, 9 


Central 


1 


17 


1 


12 


31 


0.72 


0.92 


0.42 


10 


GCR Total 


5 


50 


1 


36 


92 


0.67 


0.97 


0.40 





N ote. — Refer e nces: (1) This work; (2llMauerhan etal] J2010all; (3llCotera et al. I l ll99Sh; (4llHomeier etall j2003h ; 
(5) iMuno et all J2006h: (61 iMikles et all j2006h ; (7) iFiger et all J2002h; (81 iLiermann et all ll2009h ; (9) 
iFieer et all lll999l l; (101 jpaumard et all i2006l l; iMartins et all i2008h ; (12) iHvodo et all feOOSl) . Note that some Stars 
11, 12, and GO. 10+0.02 , shown in Figure 8, have been designated as "isolated" stars even though they are likely to 
be associated with the Arches cluster. 




2.05 2.10 2.15 2.20 2.25 

Wavelength (fim) 



Fig. 1. — f^-band spectra of confirmed 04-61 stars G359. 973— 0.008 (Star 6) and 
G359. 925— 0.049 (Star 7), which are counterparts to Chandra X-ray sources CXOGC 
J174531.4-285716 and CXOGC J174537.9-290134, respectively. 
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Fig. 2. — i^-band spectra of confirmed WNh and 04-6If + stars. 
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Fig. 3. — i^-band spectra of P Cyg-type supergiants and B supergiants. 
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Fig. 4. — i^-band spectra of confirmed strong-lined WN stars. 
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Fig. 5. — i^-band spectra of confirmed WC9 stars. 
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12 3 

H - K s 

Fig. 6. — Near-infrared color-color diagram of new WN, WC, and OB stars (green squares, 
red circles, and blue triangles, respectively). The small dots are random field stars. All of 
the massive stars exhibit infrared excess, mainly attributable to free-free and line emission 
from their dense, ionized winds. The two WC9 stars farthest from the locus of field stars, 2 
and 4, may be sources of an additional excess component of thermal emission from hot dust. 
A color version of this figure is available in the online journal. 
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Fig. 7.— Plot of observed line strength (F187N/F190N) versus i^-band SIRIUS magnitude 
for unidentified candidate Pa-excess sources (crosses), and for confirmed massive stars. The 
brightest known sources (K s < 8 mag), which are saturated in the SIRIUS catalog, were 
plotted using their 2MASS photometry. The figure includes objects from the Arches, Quin- 
tuplet, and Central Parsec clusters, and from the GCR field. Stars are marked according to 
their approximate spectral type, defined in the legend at the upper right of the figure. The 
brightness distribution of unidentified sources suggests that the sample of WN types in the 
survey might be near completion. 
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Fig. 8. — HS T /NICMOS Pa survey image from I Wang et al. 1(120 lOl ) of the Arches cluster 
region and nearby WNh stars. The marked stars 11 and 12 are new discoveries, while 
GO. 10+0.02 was previously reported in ICotera et al. I ( )1999l ). The surrounding stars all lie 



within 1-2 pc, in projection, of the cluster (beyond the tidal radius of ~ 1 pc), and they 
may have escaped from it. North is up and east is to the left in this image. 
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\ 

1 ARCMIN (-2.3 PC.) GO.120-0.048 (LBV) / 



Fig. 9.— iZST/NICMOS Pa survey image of the Sickle region from IWang et al. IfcOlOh 



which contains the Quintuplet cluster and nearby isolated massive stars, including the 06If + 
X-ra y source XI 74517 .0 from iMaue rhan. Muno. fc Morris I fl2007f ). the LBV GO.120-0.048 
from iMauerhan et al.l (l2010al ). and the new WN5b star 17 of this work. All of the unmarke d 
point sources are known massive stars spectroscopically identified in iFiger et al.l ( 119991 ) . 
North is up and east is to the left in this image. 
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Fig. 10. — HS T /NICMOS Pa survey image of the Galactic center region from lWang et al. I(l2010f ). The positions of the 
isolated massive stars are marked (i.e., those located outside of the known extent of the Arches, Quintuplet, or Central 
clusters) and color coded according to spectral type (see the image legend). The discoveries of this work (circles) are 
marked separately from the known sources (crosses). The two foreground O stars are labeled "FG". Bona fide cluster 
members are not marked. 



